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ABSTRACT 

According to one model, high-frequency quasi-periodic oscillations (QPOs) can be 
identified with inertial waves, trapped in the inner regions of accretion discs around 
black holes due to relativistic effects. In order to be detected, their amplitudes need 
to reach large enough values via some excitation mechanism. We work out in detail a 
non-linear coupling mechanism suggested by Kato, in which a global warping or eccen- 
tricity of the disc has a fundamental role. These large-scale deformations combine with 
trapped modes to generate 'intermediate' waves of negative energy that are damped as 
they approach either their corotation resonance or the inner edge of the disc, resulting 
in amplification of the trapped waves. We determine the growth rates of the inertial 
modes, as well as their dependence on the spin of the black hole and the properties 
of the disc. Our results indicate that this coupling mechanism can provide an efhcient 
excitation of trapped inertial waves, provided the global deformations reach the inner 
part of the disc with non-negligible amplitude. 
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1 INTRODUCTION 

In the past two decades, and mainly thanks to NASA's Rossi 
X-ray Timing Explorer (RXTE), about 20 X-ray binaries, 
believed to contain a black hole, h ave been analysed in detail 
iRemillard fc McClintoc^l2006al ). An importa nt character- 
istic of these objects is rapid X-ray variability (jvan der 
I2OO6I ). in which quasi-periodic oscillations (QPOs) are in- 
cluded. Of particular relevance are high-frequency quasi- 
periodic oscillations (HFQPOs), features observed in the 
power spectra of the light curves of black hole candidates, 
which are a potentially important tool in the study of strong 
gravitational fields. The frequencies observed, and their 
stability against luminosity variations, suggest a connec- 
tion w ith the inner accr etion flow l|Re miIlard & McClintock 
2006bl'). Some autho rs ("Nowak et all 11997. : .Wagoner, ,199a : 



Nowak fc Leh3 1 19991 : |kato 200 1]) argued that these oscil- 



lations can be explained in terms of modes arising in the 
accretion disc that surrounds the black hole. Turbulent vis- 
cosity, characteristic of accretion flows, generally prevents 
waves from propagating across the disc to form coherent 
global modes. A possible way for such modes to exist is for 
them to be trapped in a small region in the inner part of 
the disc. Radial trap ping of waves was first predicted by 
iKato fc Fuku3 l|l980l ). who considered oscillations in a disc 



* E-mail: B.T.Ferreira@damtp.cam.ac.uk 



around a Schwarzschild black hole, and was studied further 
by Okazaki et al... (1987. ) . Although trapped waves may re- 
veal relatively little about the properties of the accretion 
disc itself, they are a pro mising tool in the study of the cen - 
tral object. According to iRemfllard fc McCUntockl (|2006bl ). 
a good model of HFQPOs would provide the most reliable 
avenue for measuring black hole spins. Once their masses 
are also determined, an important step in testing the Kerr 
metric will be taken. 

For trapped oscillations to explain HFQPOs, an excita- 
tion mechanism for these modes is required, as their ampli- 
tudes need to reach values high enough to allow detection. 
A possibility can reside in the interaction between waves 
in t he disc and a glo bal deformation (warping or eccentric- 
ity). ICoodma^ ([l99i) initially studied the local excitation of 
waves in a tidally distorted disc via the parametric instabil- 
ity. The idea of using a warp as a n excitation m echanism 
for disc oscillations goes back to IPapaloizou fc TorquemI 
(f 19951 ). who mentioned the possibility of parametric gen- 
er ation of inertial waves. Detailed calculations are reported 
bv lGammie et all l|2000l ). More recently, and using a differ- 
ent approach, this problem was studied analytically for thin, 
relativ istic disc s with non-r otating central objects by iKatol 
(|2004 ) (see also lKatd |20o3), where similar calculations are 
made for eccentric discs). He made simple estimates for the 
growth rates of trapped inertial modes, which are of consid- 
erable interest. On the other hand, there are many uncer- 
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tainties in his calculations and he did not discuss the origin 
or nature of the global deformations. 

In this paper we develop and generalize Kato's ideas 
on this excitation mechanism and make detailed numerical 
calculations of the modes and growth rates for rotating black 
holes. We include a dynamical treatment of the warp or 
eccentricity but defer to a second paper a broader discussion 
of the origin and global propagation of these deformations. 

In Section 2 we review the trapping of inertial oscilla- 
tions in a simple, pseudo-relativistic disc model. In Section 
3 we describe the excitation mechanism for trapped inertial 
modes, which relies on a coupling between waves in the disc 
and global deformations. In Section 4 we discuss the depen- 
dence of the inertial modes' growth rates on disc parameters 
and black hole spin. Conclusions are presented in Section 5. 



2 TRAPPED INERTIAL OSCILLATIONS 
2.1 Basic equations 

The trapping of oscillations can be easily understood by 
analysi ng the fluid equations i n a s i mple isothermal disc 
model (jLubow fc Pringlelll993l : |k^[2001). Although this 
trapping happens only in discs around compact objects, a 
fully relativistic model is not necessary. The most important 
effects can be included by supplementing a Newtonian treat- 
ment with the correct relativistic expressions for the char- 
acteristic frequencies in the disc (Kato 200i|). For simplicity 
and clarity, we adopt this pseudo-relativistic approach and 
consider a strictly isothermal disc with a ratio of specific 
heats 7 = 1. Ignoring viscosity and magnetic fields, the fun- 
damental hydrodynamic equations can be written as 



imh' 



du 
'dt 



+ u-Vu = -Vh - V<&, 



-c,V ■ u. 



(1) 



(2) 



where h = log p is the enthalpy and Ca is the constant 
sound speed in the disc. We neglect self-gravitation and con- 
sider a fixed axisymmetric gravitational potential ^{r,z), 
where (r, (fi, z) are cylindrical polar coordinates. 



2.2 Linearized equations 

As in the case of stars jChristensen- Dalsgaard|[2002h . to 
study oscillations one needs to analyse what happens to this 
system of equations when the velocity and enthalpy are per- 
turbed: q — qo + q' ■ The equilibrium state of the disc is 
independent of time and azimuth: uq — ft x r — rn{r)e^, 
where Q is the angular velocity, which is independent of z 
in a strictly isothermal disc, while ho{r,z) satisfies V^o = 
rQ,^er — V"!?. Thus the perturbations acting on it can be 
written as 

q'(r, 0, z, t) — Re ^q'{r, z) exp(im0 — icijt)j , (3) 

where m is the azimuthal mode rmmber and uj is the oscil- 
lation frequency. Dropping, for simplification, the tildes and 
zeros, the linearized equations for the perturbed quantities 
can be written in the form 



-luju,. — Ziluj, = — —, 
or 



+ 2?7'"'' ^ 

dh' 
oz 

—luon — \l,zu, — — c. 



1 d{ru'j.) 'irnu'^ du'^ 
r dr r dz 



(5) 
(6) 
(7) 



where lu = uj — mO, is the Doppler-shifted wave fre- 
quency, which is zero at the corotation radius, and k and 
flz are the epicyclic and vertical frequencies, respectively 
l|Binnev fc TremainelllQSSi ). We apply the thin disc approx- 
imation {dh/dz = — n^z), and neglect the term u',.dh/dr in 
the last equation. The latter approximation is valid if the 
radial wavelength for perturbations is smaller than the ra- 
dial scale on which the enthalpy varies in the basic state. 
Variables can then be further separated in r and z, using 
(Okazaki ct al. 1987) 



{u'^,u'^,h') = {ur{r),u^{r),h{r))Re„ (^^^ 



M2(r)IIe„_ 



(8) 
(9) 



where Hen is the modified H ermite polynomial of order n 
(jAbramowitz fc Stegunlll972l '). with n = 0, 1,2 , 3, . . . being 
the vertical mode number and H — ^ c^jflz the vertical 
scaleheight of the disc. (Since He_i is not defined, for n = 
we have u.^ =0.) It should be noted that this separation of 
variables is not exact since H depends on r (as described by 
JKato (200l|), this separation is valid to lowest WKB order; 
[Nowak fc Wagoner] jl992h use a slowly varying function of 
r to separate variables). The variation of H with r couples 
different vertical modes l|Tanaka et al.|[2002l ) but this effect 
is weak when the radial wavelength is short, and we neglect 
it. The final set of ordinary differential equations in r for the 
perturbed quantities reads 



— 2VLu^ — — — 



+ 2?i"'- 



dft 
dr ' 
imh 



-iuih — fl.Huz ~ 



2.3 Wave modes 



1 d{rur) imu^ 



dr 



+ 



(10) 

(11) 
(12) 

(13) 



The dispersion relation for wave modes in the disc can be 
determined by further assuming that the radial wavelength 
of the perturbed quantities is much smaller than both the 
azimuthal wavelength and the characteristic scale for radial 
variations of the equilibrium quantities. It can then be veri- 
fied that perturbatio ns with local radial wavenumber k obey 
(jOkazaki et al.lll987r ) 



(14) 



As argued before, the most important relativistic effects on 
wave propagation can be included by using relativistic ex- 
pressions for the characteris tic frequencies. For a particle 
orbit, these read (|Katolll990l ). 



(4) 



Q = (r=*/' + a) 



(15) 
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n /i 6~ 8a 3a2 
'' = "V r + ;^"~' ^^^^ 

^ ^ 4a 3a2 

"^ = "V^"^ + ~' ^^^^ 

where a is the dimensionless spin parameter of the central 
object (—1 < a < 1), r is in units of the gravitational radius 
Tg = GM/c? and the frequencies are in units of c^/GM, 
where c is the speed of light, G the gravitational constant 
and M the mass of the black hole. We assume that the 
characteristic frequencies in the disc can be approximated 
by these particle orbit expressions. Of great importance is 
the variation of the epicyclic frequency k with r: as r de- 
creases, K increases, reaches a maximum and then goes to 
zero at the radius of the marginally stable orbit, rms, which 
can be regarded as the inner edge of the accretion disc. Since 
we are interested in studying oscillations in accretion discs 
around compact objects, we adopt these relativistic expres- 
sions. They correctly describe the frequency and stability of 
orbits in the Kerr metric, as well as the apsidal and nodal 
precession rates, but some information about the metric co- 
efficients is lost in this pseudo-relativistic approach. 

From the dispersion relation we can see that, if n 7^ 
0, two types of wave-like solutions, propagating in different 
regions in the disc, are possible: a high-frequency one with 
Cj^ > raa.y^{K? ,nQl) = nfi^ (p mode) and a low- frequency 
one with uj'^ < mm{K^ , nfil) — (r mode). If n = the 
dispersion relation for the non-trivial mode becomes uj^ = 
k^Cs+K-^ and waves can propagate where li^ > k^; this is the 
inertial-acoustic mode. In non-isothermal discs the closest 
equivalent of this mode behaves like a surface gravity wave 
or steUar f mode (|Ogilvielll99i iLubow fc Ogilvielll998l ) but 
here we refer to it as the n = mode or 2D mode, since it 
involves a purely horizontal motion independent of z. 

The inertial o r r modes (so called by 

iKorvcanskv fc Pringld Il995t ) are nearly incompressible 
since they are restored by inertial forces, avoiding acoustic 
effects. (They are often called g modes in the literature, 
but they are not related to internal gravity waves or stellar 
g modes.) Acoustic or p modes have pressure as their 
main restoring force and are essentially compressible. In 
an isothermal disc, the 2D mode is a purely horizontal 
compressible mode. 

The propagation regions of the p and r modes allow us 
to define important rad ii in the disc, the resonant radii (e.g. 
iLubow fc Qgilvielll998l ). Non-axisymmetric waves can have 
three types of resonances: corotation where li = 0, Lindblad 
resonances where — k^ — 0, and vertical resonances where 
(jj^ — nQ^ = (for n 7^ 0). Lindblad resonances are turning 
points for the r and 2D modes, while vertical resonances are 
turning points for p modes. 

The three different types of modes propagate in different 
regions in the disc. If the epicyclic frequency has a maximum 
at some particular radius, r and 2D modes can be trapped in 
the inner part of the disc, while p modes always propagate 
to the out er boundary beyond t he outer vertical resonance. 
Originally iKato fc Fukud (|l980l ) considered the trapping of 
2D modes in the very inner region of the disc. These modes 
can be trapped between the radius of the marginally stable 
orbit and the inner Lindblad resonance. However, the condi- 
tions at are not well understood and it is not clear if this 
trapping region can work as a resonant cavity. Therefore, 



here we focus on the trapping of r modes, which happens 
below the maximum of the epicyclic frequency between two 
Lindblad resonances (and with no corotation resonance in 
between), a resonant cavity naturally created by the non - 
monotonic variation of k with radius (|Okazaki et al.|[l987l ). 

Of particular importance is the axisymmetric trapped 
wave with frequency uj « max(K), and with the simplest 
possible radial structure. This mode is trapped in a small 
region close to the maximum of the epicyclic frequency, and 
is important, not only because it is naturally confined and 
therefore more likely to occur in the presence of turbulent 
viscosity, but also because its frequency can be identified 
with max(«:), which depends only on the properties of the 
black hole: its mass M and its angular momentum a. Fur- 
thermore, this mode is likely to be most easily observed, 
as it may produce a net luminosity variation of the disc 
without cancellations. Therefore, measuring the frequency 
of this mode and determining the mass of the central object 
by, e.g., studying the orbit of its binary companion, one can, 
in principle, find the spin of the black hole. 



2.4 Numerical calculation of trapped r modes 

To find the radial structure of trapped modes we need 
to solve the system of equations pup - (ll3|l . subject to ap- 
propriate boundary conditions. Numerical calculations of 
waves trapped near the maximum of th e epicy clic frequency 
were first performed bv lOkazaki et all (|l987l ). Here we fo- 
cus on the simplest possible trapped inertial modes, with 
m = and n — 1. An approximate analysis of equa- 
tions (|10p - (|13|) close to the maximum of the epicyclic fre- 
quency shows that, between the two Lindblad resonances, 
these solutions are described by parabolic cylinder func- 
tions (jAbramowitz & Stcgun 1972) involving Hermite poly- 
nomials of order I — 0, 1, 2, ... , centred at the maximum 
of K , like the solution s of the quantum harmonic oscilla- 
tor (|Perez et al.l Il997l fl. The lowest order mode {I = 0) 
has a Gaussian structure in r. In this section we solve the 
same problem using numerical methods and fewer approxi- 
mations, as a prelude to an analysis of the non-linear mode 
couplings that cause these modes to grow. 

Since we expect to find trapped modes only for some 
discrete values of the oscillation frequency, we solve the set 
of equations as a generalized eigenvalue problem, of the form 

AU = -iioBU, (18) 

where U is the column vector whose components are the 
r mode quantities {ur,u^,Uz,h) evaluated at a set of dis- 
crete points, A is the matrix representing the system of 
equations (|10p - (|13[) . and B can be different from the identity 
matrix depending on the boundary conditions used. To solve 
this problem numerically, we use a pseudo-spectral method 
with Chebyshev polynomials. We use a Gauss-Lobatto grid, 
x{i) = cos{ni/N), where A*" is the number of grid points, and 



^ It should be noted that iPerez et al. I lll997h use n as the ra- 
dial mode number. Here the three quantum numbers are (Z,m,n) 
corresponding to the three coordinates (r, </), z). 
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the Chebyshev coordinate 
dial coordinate by 



-1 < X < 1 is related to the ra- 



gout 



+ 



(19) 



so that Tin < r- < rout, where r in — fms and Tout IS an outer 
radius chosen to be larger than the outer Lindblad resonance 
for the r mode. The representation of the first derivatives in 
the matrix A is achieved usi ng the Cheb yshev collocation 
derivative matrix, defined by l|Bovdll200ll) 



'(l + 2iV2)/6 
-(1 + 27V2)/6 



i=j = 
i=j = N 
i = j, j^O,N 



(20) 



where po ~ Pn ~ 2, and pj — 1 otherwise. 
The generalized eigenvalues and eigenfunctions of A are 
then calculated numerically, using IDL's eigenvalue solver, 
LAJEIGENPROBLEM, which uses a QR decomposition, 
and is based on LAPACK routines. The boundary condi- 
tions used are the following: 

At rin Ur = 0. According to the dispersion relation, 
the r mode is expected to be exponentially decaying for radii 
smaller than its innermost Lindblad resonance and there- 
fore its velocity is supposed to be approximately zero at 
the marginally stable orbit, which justifies the choice of this 
boundary condition. 

At rout dur/dr — ikur, where k is given by the dis- 
persion relation (|14|l for m = 0, n = 1 at rout, and using 
Lu ~ max(K). Since the r mode can propagate again as a 
p mode outside the vertical resonance, we choose an outer 
radius beyond this location, so that the mode is oscillatory 
there, i.e., k is real, and we select the outgoing wave solution 
by choosing fc > so that the group velocity of the waves 
at rout is positive. This condition allows the wave to lose 
energy through the outer boundary and minimizes artificial 
wave reflection there. 

We allow for the possibility that uj is complex, in which 
case its imaginary part is the growth rate of the disturbance. 
In fact, in the absence of non-linear mode couplings, we ob- 
tain slowly decaying solutions with Im(aj) < as a result of 
the outgoing-wave outer boundary condition. 

In Fig. [T] we show the variation of Ur with radius for 
two typical trapped solutions, corresponding to two differ- 
ent radial mode numbers I = and / = 1. (Since we are 
solving an eigenvalue problem, solutions are multiplied by 
an arbitrary amplitude.) The complex frequencies of the 
modes represented in Fig. [T] are 0.03196 - 3.6884 x 10"^i 
and 0.02989 - 1.1846 x lO^'^i (in units of c^/GM), respec- 
tively. In these units, and for the value of a used, the maxi- 
mum of K is 0.03312. These modes are slightly damped be- 
cause of the boundary condition used at the outer radius, 
which selects the outgoing wave only. If the sound speed is 
smaller, the frequency of the modes is closer to the value of 
the epicyclic frequency at its maximum, and the trapping 
region and damping rate are smaller (see Tab. [T|. We ob- 
tained identical results by solving the two-point boundary- 
value problem using a shooting method. 

These inertial modes can be thought of as waves 
trapped in a virtual potential, U(r) = ~k{r)^ (|Li et al.l 



Cs/C 


frequency 




0.002 


0.03289 + O.Oi 




0.005 


0.03254 - 1.276 x 10" 


10; 


0.01 


0.03196 - 3.688 x 10 




0.02 


0.03079 - 2.417 x 10 


~5i 



Table 1. Dependence of the real and imaginary part of the Z = 
r mode frequency (in units of c^/GM) on the sound speed Cg, for 
a = 0.5. We use 150 collocation points and a value of 18.2331 for 
the outer radius, in units of GMjf? (rjn = rms = 4.2331). 
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Figure 1. Variation of the real part of the radial component of the 
axisymmetric, rt = 1 r mode velocity with radius for Cs/c = 0.01 
and a = 0.5 for (a) I = and (b) 1 = 1. The triangle indicates the 
radius where the epicyclic frequency is maximum, and the crosses 
and asterisks the Lindblad and vertical resonances, respectively. 



l2003h . which for a frequency close to the maximum of k 
is similar to the harmonic oscillator potential. If U{r) < 
waves can propagate, being evanescent in the regions where 
potential barriers exist. Also, as in quantum mechanics, 
these trapped inertial waves can escape through the poten- 
tial barriers and propagate on the other side, as p modes. 
The inertial modes are evanescent between the inner ra- 
dius and the flrst Lindblad resonance and between the sec- 
ond Lindblad resonance and the vertical one. The 'leak- 
age' through the potential barrier is also verifled, as our 
results show small-amplitude oscillations after the vertical 
resonance (Fig. [l]). The larger the sound speed, the larger 
the width and smaller the height of the barrier, and more 
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'leakage' through the barrier is expected, which is verified 
numericaUy (Tab, [l] ). This 'leakage' was first predicted by 
lOkazaki et all (jigSTf ). 

As described above, the inertial mode characterized by 
(I, m, n) — (0, 0, 1) is likely to be relevant to the interpreta- 
tion of observed oscillations. In the next section we describe 
an excitation mechanism for this mode, based on non-linear 
wave coupling. If the modes are not strongly coupled, we ex- 
pect the mechanism to affect mainly the growth rate Im(aj) , 
so that the structure of the wave still resembles a Gaussian 
centred at the maximum of k, as shown in Fig. [T](a). 



3 GROWTH OF OSCILLATIONS IN 
DEFORMED DISCS 

Low-frequency modes with azimuthal mode number m = 1 
have been widely stu died in the context of (quasi-)Keplerian 
accretion disc theory (iKatdfToSS ) . since they are global, typ- 
ically vary on length-scales c omparable t o the radial extent 
of the disc, and are long-lived. iKatol l| 19891 ) showed that these 
modes are also global in a relativistic disc. A global m = 1 
mode with one node in the vertical direct ion (n = 1) is typ- 
ically identified with a warp in the disc (|Papaloizou fc LinI 
1 19951 ). while n — modes correspond to eccentric discs. This 
is easily seen if we focus on the action of each of these modes 
on a ring. The vertical displacement of a m = 1, n = 1 mode 
is independent of z, and proportional to cos{(j} — constant) 
at fixed r and t, which corresponds to a tilting, as the dis- 
placement with respect to the disc plane is different at each 
azimuthal angle. In the case of the n — mode, the radial 
displacement is the one that is independent of z, and propor- 
tional to cos{(j} — constant) at fixed r and t. This means that 
the displacement within the disc plane varies with (f>, cre- 
ating an elliptical orbit. As we discuss below, these global 
deformation modes can be produced by the presence of a 
binary companion or by instabilities. 

3.1 Warped discs 

It is believed that warps exist in many astrophysical discs. 
Precessing warped discs have been used successfully to ex- 
plain long-term light- curve variations in Her X-1 and sorrt e 
other X-ray binaries (|Kat j|l973l : iGerend fc Bovntonlll97^ . 
If the central object is a compact radiation source, warps 
can be induced by rad iation pressure forces (jPringlelll993 : 
IWiiers fc Pringld il999l ). However, this mechanism is less 
likely to operate in systems containing a b lack-hole primary 
because the disc needs to be very large (jOgilvie fc Dubua 
I2OOII ). 

If the central object is a rotating black hole, its axis 
of rotation might not be perpendicular to the plane of the 
binary orbit in which the accretion disc forms, in which 
case the disc is said to be misaligned or tilted. There is 
both th eoretical and observ ational evidence for this tilt- 
ing (see [Fragile et al.l |20o3) and references therein). The 
misalignment of the orbital angular momentum of the disc 
and the spin angular momentum of the black hole re- 
sults in important changes in the structure of the in- 
ner disc as it will be subject to Lense-Thirring preces- 
sion (jBardeen fc PettersonI Il975l ) . This differential preces- 
sion tends to twist the disc, which may adopt a station- 



ary warped shape. Depending on the 'amount of viscos- 
ity', the induced war ps can propagate ei t her d iffusively, 
roughly speaking if the lShakura fc SunvaevI () 19731 ) viscosity 
parameter a is greater than H /r, or in a wave-lik e man - 
ner if a is smaller than H/r. Ilvanov fc Illarionovl (|l997l ) 
showed that the warp has an oscillatory radial structure in 
a low- viscosi t y disc , and this was investigated further by 
iLubow et ai] (|2002l ). 



3.1.1 Variation of disc tilt with radius 

In a vertically isothermal (pseudo-)relativistic disc, with 7 — 
1, a zero-frequency mode with n = 1, m = 1 can propagate 
at all radii if a > (i.e. if the disc and black hole rotate 
in the same sense). This can be seen from the dispersion 
relation (|14p and the expressions (|16p and (|17|l . because in 
this case lo^ — is greater than both and Q,^, and 
therefore k > 0. Such a stationary warp can be described 
by linear perturbations of the form 



{u'wr,u'v^^,h'w) = {uwr{r),uyf4,{r),hw{r)) z, 

ttWz = WWz(r), 



(21) 
(22) 



where the subscript W refers to warp quantities, and the 
dependence e"^ is understood. The simplest possible warp 
solution is the rigid tilt, valid for a non-rotating black hole 
{a — 0, Q = Qz), described by uwz = WQr, uwr = —Wil, 
uw4> ~ ~'iWd{rQ,)/dr and h w = —iWQ^r, whe r e W i s the 
constant tilt inclination (see iPapaloizou fc LinI (|l995l ) but 
note that they use g instead of W to represent the disc tilt). 
If W varies with r, equations ©-([7]) with m = 1, lu — 
admit a solution of the form 



UWr 



-Q.Wz 4- rz 



AW Q.^ 



,^^ = _iW^,_(,n)+i^^ 

/iW = -if^^Wrz, 

where Wir) is the solution of 



m 

dr ' 



d_ 

dr 



dW 



dr 



1 m 

r dr 



-w. 



(23) 

(24) 

(25) 
(26) 

(27) 



This equatio rEI is closely r e lated, but not identical, to equa- 
tion (17) of iLubow etal] (|2002l ). which was derived from 
an analysis of global warps in discs that are not necessar- 
ily isothermal. We solve this equation numerically, using a 
4th order Runge-Kutta method with the boundary condi- 
tion dM^/dr(rin) = 0, corresponding to zero torque at the 
inner edge. The amplitude of this linear solution may be 

^ The above analysis uses the relation = 4f2^ -|- 2rf2df2/dr, 
which is not exactly true of the relativistic expressions because 
r^O is not quite the specific angular momentum in relativity. On 
the other hand, since the pseudo-relativistic treatment is not fully 
self-consistent, if this Newtonian relation is not used here, the 
rigid tilt solution is not obtained for a = 0, contrary to what is 
expected physically. Therefore, we choose to use the Newtonian 
relation in our treatment of the warp. This is not expected to 
significantly influence the final results. 
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Figure 2. Warp function W{r) for a = 0.5. The sound speed is 
0.01c and Winn) = Wo = 0.003. 



fixed by specifying tiie value Wo = W(rin) at the inner 
boundary, i.e., at tlie marginally stable orbit; this corre- 
sponds to the (small) inclination of the inner edge of the 
disc with respect to the equator of the black hole. A typ- 
ical solution is shown in Fig. [2l The warp ha s an o scilla- 
tory behaviour, as found bv llvanov fc Illarionovll|l99'i1 ). with 
the wavelength increasing with radius, consistent with the 
local dispersion relation. This non-monotonic behaviour of 
the inclination contrasts wi th the Bardeen-Petterson effect 
(|Bardeen fc Petterson|[l975l ) , which was derived using an in- 
correct equation for the warp. We would normally expect 
W{r) to tend to a constant value at large r, corresponding 
to the inclination of the outer part of the disc with respect 
to the equator of the black hole. Unfortunately this is not 
true of the approximate equation (|27p . which does not hold 
accurately at large r because the wavelength becomes com- 
parable to the radius. However, since we are interested in the 
interaction of the warp with waves that propagate in the in- 
ner disc, this is not expected to significantly affect the final 
results. We defer to a second paper a more realistic treat- 
ment of the propagation of the warp into the inner part of 
the disc. 



3.1.2 Coupling mechanism 

The non-linearities in the basic equations ([T| and ((5} provide 
couplings between the different linear modes of the system. 
We are interested in those couplings that lead to amplifica- 
tion of the trapped modes. T he basic ide a of the excitation 
mechanism in warped discs (|Katol [20041 ) is that the warp 
interacts with a wave in the disc (the trapped r mode) giv- 
ing rise to an intermediate mode. This intermediate mode 
can then couple with the warp to feed back on the original 
oscillations (see Fig. [3]), resulting in growth of the latter. 

For the r mode to be excited, it needs to gain energy in 
this coupling. Since the warp has null frequency, its energy 
is essentially zero and so the energy exchanges only happen 
between the r and the intermediate modes and the disc. It 
is widely agreed, and certainly true in the short-wavelength 
limit, although a general proof is lacking, that a mode that 
propagates inside its corotation radius has negative energy, 
i.e., the total energy of the disc is reduced in the presence of 
the wave, which is possible because the disc is rotating. On 



r— mode 
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Figure 3. Diagram representing the interactions involved in the 
coupling mechanism involving a warped disc 



the other hand an axisymmetric wave, such as the r mode, or 
one that propagates outside its corotation radius, has posi- 
tive energy. Suppose that, through coupling with the warp, 
the r mode generates an intermediate wave that propagates 
inside its corotation radius and therefore has negative en- 
ergy. In the process of generating this wave, the r mode 
gains energy and is amplified. For sustained growth of the 
r mode, the intermediate wave must be damped so that its 
negative energy is continually replenished by the r mode. 
(The damping process itself draws positive energy from the 
rotation of the disc.) Therefore a dissipation term should 
be included in the equations for the intermediate mode. We 
choose to damp this wave locally at a rate f3i}, where /3 is 
a dimensionless parameter. The origin of this term is not 
discussed here but if we interpret it as some type of viscous 
dissipation or friction in the disc we expect the intermedi- 
ate mode, which propagates in a larger region in the disc, 
to be more affected by it than the r mode, as the latter is 
trapped in a small region, and has a simpler radial struc- 
ture. Also, the intermediate mode approaches its corotation 
radius (or the marginally stable orbit), where it is expected 
to be absorbed, and this effect is implicitly included in the 
intermediate mode equations when the friction term is in- 
cluded. Therefore, we neglect the dissipation term in the 
equations for the r mode. The growth rate that we obtain 
for the trapped mode should be compared with estimates of 
its damping rate due to turbulent viscosity. 

For the coupling to occur the waves need to propagate 
in the same region in the disc and the parameters u; and m 
for the 3 modes need to follow some basic coupling rules. 



(28) 



where the subscripts R, W and I refer to r mode, warp 
and intermediate mode quantities, respectively. These rules 
follow from the quadratic nature of the non-linearitites in 
the basic equations ((TJ and Also, we can get informa- 
tion about the vertical mode number m of the interme- 
diate mode by remembering that the vertical dependence 
is given by Hermite polynomials. The warp quantities are 
proportional to z and the simplest possible r mode has 
one node in the vertical direction, therefore its quantities 



On an excitation mechanism for trapped inertial waves 7 



are proportional to Hei ~ z. When the warp and this 
r mode interact, the couphng terms will be proportional to 
/H^ = {z^ IH^ - 1) + 1 ~ He2 + Heo, i.e., they give rise 
to two intermediate modes, one with 2 nodes in the vertical 
direction, rii — 2, and a 2D mode with m = 0. We consider 
both possibilities. Note that the frequency of these inter- 
mediate waves is that of the r mode, and they are present 
because they are forced through the couplings. They could 
not exist as free waves satisfying the boundary conditions 
at this frequency. 

According to rules (I28p . if ujr — uj, then loi = uj. For 
the azimuthal mode numbers, we have mi = ttir ± 1. By 
analysing the propagation regions for these modes, and con- 
sidering u! ~ max(K + mil), it is possible to conclude that 
the n — mode with mi — mR — 1 and frequency lu does 
not have any Lindblad resonance, i.e., the point where the 
wave should be excited, in the disc. On the other hand, the 
mode with mn + 1 and frequency uj has an inner Lindblad 
resonance close to the region of propagation of the r mode. 
For these reasons, we choose the azimuthal mode number 
of the intermediate mode to be ma -I- 1 = 1, if the r mode 
is axisymmetric. This mode propagates inside its corotation 
resonance, while the r mode has positive energy. The pres- 
ence of an inner Lindblad resonance means both that the 
intermediate mode attains a larger amplitude than it would 
in the case of nonresonant forcing, and that the flow of en- 
ergy is such as to amplify the r mode. 

If the vertical mode number of the intermediate mode 
is chosen to be 2 instead of 0, the r mode with {uj,m,n) = 
(a;,0, 1) interacts with a (lj,1,2) intermediate mode. The 
former propagates where lo'^ < , while the latter propa- 
gates where {ui — Q,)^ < , where ui is slightly less than 
max (k). The propagation regions overlap close to the max- 
imum of the epicyclic frequency since « 2k in this region. 

The interaction of the r mode with the n = 2 inter- 
mediate mode in a warped disc must be treated carefully 
because in this case the intermediate mode propagates be- 
tween its Lindblad resonances and is absorbed at the coro- 
tation resonance. This is a radius in the disc where the po- 
tential U{r) — —k^ tends to infinity and is therefore difficult 
to treat numerically because the wavelength tends to zero. 
One way of solving this problem is by including a relatively 
strong dissipation term in the equations for the intermediate 
mode. In this way the wave excited at the inner Lindblad 
resonance is damped before reaching corotation. This also 
works for the energy exchanges between the modes and the 
disc, since the n = 2 mode has negative energy in the region 
where it is damped. 

3.1.3 Results 

To find the r mode growth rate resulting from these interac- 
tion we solve the systems of coupled equations (|A19|I - (|A25|) 
and (|A26|) - (|A33|) (see Appendix A) for the interactions of 
the r mode with the n = and n = 2 intermediate mode 
respectively. By considering the warp to have a fixed ampli- 
tude and neglecting the feedback of the r mode and inter- 
mediate modes on the warp, we still obtain a linear system 
of equations, although now the r mode and intermediate 
modes are coupled through the warp. We treat the n = 
and n = 2 intermediate modes separately, although in prac- 
tice both coupling mechanisms act simultaneously and the 
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Figure 5. Variation of the real part of the radial component of 
the (a) m = 1, n = 0, and (b) m = 1, n = 2 intermediate mode 
velocity with radius for Cg/c = 0.01, a = 0.5, W{ri^) = 0.003, 
and P = 0.1. The triangle indicates the radius of the corotation 
resonance, and the crosses the Lindblad resonances. 



net growth rate is the sum of the rates due to the individual 
mechanisms. 

We solve these systems numerically, using the Cheby- 
shev method described in Section^ For the r mode the same 
boundary conditions as before are used. Similar conditions 
are applied to the intermediate modes, i.e., uir = at rin 
and Auir/Ar = ikiuir at rout, where fci is given by the dis- 
persion relation (|14|) at rout for m = 1 and n = or n = 2, 
depending on the intermediate mode we are considering. As 
for the r mode, we choose the sign of ki so that the outgoing 
or exponentially decaying wave at rout is chosen. The choice 
of the inner boundary condition for the radial component 
of the velocity of the n = 2 mode is justified by the fact 
that this mode is exponentially decaying there. This choice 
is harder to justify for the n = mode since it is oscillatory 
at rin. This means that if uir — there then the wave is 
reflected at the inner boundary. Since the conditions at the 
marginally stable orbit are not clear, we cannot be sure of 
the physical validity of this condition; we choose it because of 
its simplicity. Rigorously we would need more boundary con- 
ditions to solve this problem, since more than 4 derivatives 
appear in each system. Ifowever, since the coupling terms 
are expected to be small, mr^ is roughly proportional to 
URr (and similarly for other quantities), thus the boundary 
conditions imposed for the latter will be indirectly imposed 
to the former. 
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The aim is to find the frequencies lo for which the so- 
lutions corresponding to the r mode are trapped, i.e., for 
which URr resembles the parabolic cylinder functions as in 
Fig. [T] (a) , which is expected if the coupling terms are small 
when compared to the other terms in the equations. The 
imaginary part of cu then gives the growth rate (or damp- 
ing rate, if it's negative) of the trapped r mode. In Fig. \^ 
we show the n = and n = 2 intermediate modes involved 
in the coupling process, when the dissipation is strong. It 
should be noted that when the dissipation is weak, the n = 2 
intermediate mode develops a very short wavelength as it 
approaches the corotation resonance. If l3 is too small, the 
length-scale on which this wave dissipates is not resolved 
by our numerical method. The variation of the growth rate 
with several parameters is shown in Fig. ID These results are 
discussed in Section 4. 

3.2 Eccentric discs 

Another possible mechanism for the excitation of trapped 
waves is their interaction with an eccentric disc ((Katai2007l ) . 
Interacting binary stars with mass ratio q < 0.3 are believed 
to have eccentric accretion discs. This phenomenon is well 
documented in the case of cataclysmic variable stars, where 
superhumps are observed during the superoutbursts of the 



SU UMa class of dwarf novae l|Patterson et al.ll2005l ) and in 
other systems of low mass ratio. In these systems, a resonant 
interaction of the orbiting gas with the tidal pot ential of the 
companion star allows a growth of eccentricity (|Whitehurstl 
11988; Lubow 1991a , b). Superhumps are also observed in an 
increasing number of low-mass X-ray binaries, and systems 
exhibiting black-hole HFQPOs are likely to have mass ratios 
q < 0.3 and therefore to have eccentric discs during at least 
some phases of thei r outb ursts. 

Recently, iKatol (|2007l ) argued that one-armed global os- 
cillations, symmetric with respect to the z = plane, can 
excite trapped oscillations. His conclusions are based on an- 
alytical, Lagrangian calculations and are too crude to allow 
for more than simple estimates for the growth rates. In this 
section we describe an excitation mechanism similar to the 
one reported previously, but where an (m = l,n = 0) ec- 
centric mode has the role that previously belonged to the 
(m = l,n = 1) warp wave. Using the same numerical 
method as before, we calculate the trapped r mode growth 
rates. 

3.2.1 Variation of eccentricity with radius 

As before, consider the set of equations (|10p - (|13|) . A global 
eccentric mode corresponds to a zero-frequency wave with 
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m = 1 and n = Q. (If the global eccentric mode precesses 
freely, the frequency is not exactly zero but is completely 
negligible compared to the characteristic frequencies in the 
inner part of the disc.) In this case, equations (|10|) " (|13|l are, 
after the usual separation of variables, reduced to 



lQ,U-Er — 2Q.Ue^ — — , 



initE</) + 



20. 
ifJuEz = 0, 

iS7/iE = — c? 



Id .me0 

--;-(rUEr) + 1 

r ar r 



(29) 

(30) 
(31) 
(32) 



where the subscript E refers to eccentric mode quantities. 
This system admits a solution of the form 



UEr = iEQ,r, 

_ dr^Q. AE 
Q?r'^ — Cg dr 

^ft dE 



-rE. 



dr 



(33) 
(34) 
(35) 
(36) 



where E{r) is the eccentricity of the disc at radius r, and 
satisfies 



1 d 



(^^-r!")i5 = -— , — 



r^c^n^ dE 



dr 



(37) 



Again, this equati on is closely related, but ri ot identical, 
to equation (21) of iGoodchild fc Ogilvid (|2006h . which was 
derived from an analysis of global eccentricity in a two- 
dimensional disc. Radially propagating solutions are ob- 
tained because < Q,^ in a relativistic disc. As for the 
warp tilt W{r), we solve this equation numerically, using a 
4th order Runge-Kutta method with boundary conditions 
E{rin) = Eo, and dE/dr{rin) = 0, where Eq is an arbitrary 
value for the eccentricity at the inner boundary, i.e., at the 
marginally stable orbit. A typical solution for E{r) is shown 
in Fig.[6l The eccentricity has an oscillatory behaviour, with 
the wavelength decreasing with radius, consistent with the 
local dispersion relation. Again, we defer to a second paper 
a more realistic treatment of the propagation of the eccen- 
tricity into the inner part of the disc. 

3.2.2 Coupling mechanism 

The excitation mechanism is similar to the one discussed 
in the previous section: a global deformation mode, which 
is now the eccentricity mode, characterized by (aj,m,n) = 
(0, 1, 0), interacts with a trapped r mode {uj, 0, 1) giving rise 
to an intermediate mode. The latter then couples with the 
global mode to feedback on to the trapped r mode. 

Coupling rules require the intermediate mode to have 
the same frequency as the trapped wave and, as before, we 
have mi = 1, in the case where the trapped r mode is ax- 
isymmetric. As for the vertical dependence, since the eccen- 
tric mode has n = 0, the intermediate mode can only have 
the same vertical mode number as the trapped mode, i.e., 
ni = 1. The propagation region for this mode is the same as 
for the (cj,l,2) intermediate mode, present in the interac- 
tion of the r mode with the warp. As in that case, a relatively 



0.004 



0.002 



^ 0.000 



u -0.002 



-0.004 




10 20 30 40 

radius in units of GM/c'2 



50 



Figure 6. Eccentricity function E{r) for a = 0.5. The sound 
speed is 0.01c and E(rin) = Eq = 0.003. 



large damping term must be included in the equations for 
the intermediate mode so that it dissipates on a resolved 
scale before reaching the corotation resonance. In this case 
the energy exchanges are similar to the ones discussed for 
the interaction in a warped disc. 

To find the growth rates that result from this interac- 
tion, we solve equations (|A34|l - (|A41|l . using the same nu- 
merical method and boundary conditions as before. 

The variation of the growth rate with the inner eccen- 
tricity, sound speed, spin of black hole and dissipation factor 
is shown in Fig. [T] The values of the growth rate achieved in 
the interaction of the trapped wave with the eccentric mode 
are, in general, similar to the values obtained in the interac- 
tion with the warp, if the inner inclination and eccentricity 
are similar. 



4 DISCUSSION 

In this section we discuss the results shown in Figs |4] and [71 
where the dependence of the I — r mode growth rate with 
several parameters is represented. 



4.1 Growth rates in warped discs 

In the variation of the growth rate with the dissipation 
factor, two regimes can be considered: weak dissipation 
iP ^ 0.05), where the variation is approximately linear, and 
strong dissipation (/3 > 0.05), where the growth rate re- 
mains approximately constant when /3 varies (Fig. [4] (a)). In 
the former case, the n = intermediate mode is launched at 
its inner Lindblad resonance and propagates, being slightly 
attentuated, until it reaches the inner boundary where it is 
refiected. Owing to the attenuation, the refiected wave am- 
plitude is smaller than the incident one, and therefore the 
intermediate mode does not cancel itself, leaving a small 
amount of energy available for the r mode to be excited. 
In the strong dissipation regime, the n — mode is dissi- 
pated before reaching the marginally stable orbit. In this 
case, all the energy carried by this wave becomes available 
to excite the trapped mode. As for the n = 2 intermediate 
mode, in the strong dissipation regime, the wave is com- 
pletely dissipated before reaching the corotation resonance. 



10 Barbara T. Ferreira and Gordon I. Ogilvie 



5x10"' 

° 4x10"= 
o 

a, 3x10-= 
"o 

5 2x10"= 
o 

x; 

I 1x10"= 



cs = 0.01, = 0.5, e(rin) = 0.003 



= 0.5, beto=0.1, cs = 0.01 



0.00 












* J 






















: 








(o) : 



0.05 0.10 0.15 

domping foctor beto 



0.20 




0.01 0.02 0.03 

eccentricity at inner rodius 



0.04 



2.0x10- 



^ 1.5x10"* 



1.0x10-* - 



I 5.0x10" 



= 0.5, beta = 0.1, e(rin) = 0.003 





3.6E-9/<-2 




















"IT - - - - - 



0.004 0.006 0.008 0.010 0.012 0.014 0.016 
sound speed in units of c 



5x10'= 



4x10 



3x10 



0) 2x10 



1x10 



08=0.01, beta=0.1, e(rin)=0.003 




0.4 0.6 
spin of BH 
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Physically we would expect no dissipation term to be neces- 
sary in this case. An arbitrarily small amount of dissipation 
should lead in principle to the complete absorption of the 
wave at the corotation resonance. However, this cannot be 
verified numerically because of the difficulty in resolving the 
wavelength of the intermediate mode as the singularity at 
the corotation resonance is approached. 

For small warp amplitudes, i.e., before the coupling 
terms start affecting the structure of the eigenfunctions, the 
growth rate grows with the square of the warp amplitude at 
the inner boundary (Fig. U (b)). This is expected since the 
coupling mechanism relies on the 'use' of the warp twice: 
first on the interaction with the r mode to give rise to the 
intermediate modes, and then again on the interaction with 
the latter to feed back on the former (Fig. [Sjl . 

The excitation mechanism discussed here is similar to 
the well known parametric instability, in the case where one 
of the modes is strongly damped. The parametric instability 
is a type of resonant coupling between three modes satisfy- 
ing (jjp ~ ujdi + ij-id2 , where the subscripts p and d refer to 
parent and daughter modes, respectively. The parametric 
instability results in the transfer of energy from the former 
to the latter, when the daughter modes have small ampli- 
tude. The equati ons describing the evo lution of the mode 
amplitudes read l|Wu fc Goldreichll200ll . adapted from), 



d^d 



dt 

dAi2 
dt 



= -7di^di — iciJdiAdi + iu)di(^ApA'^2, 

= ~7d2^d2 - il^d2^d2 + ilJJd2(^A*iiAp, 



(39) 
(40) 



where 7j > is the linear amplitude growth/damping rate 
of mode j and a is the non- linear coupling constant. Let 
us consider a simplified case where the amplitude of the 
parent mode is approximately constant in time (because the 
daughter modes are of small amplitude), and Udi = — 'i'd2 = 
7p — 7di = and 7d2 ~ 7- In this case, the parent mode 
can be compared to the warp while the daughter modes 1 
and 2 can be compared with the r and intermediate modes, 
respectively. Assuming A^i oc exp (st), the growth rate is 



Re(.) = -|+ [^ + |Ap|V^^' 



(41) 



(38) 



If 7 <S IApIcto;, Re(s) ~ lAplao; — ^, i.e., the growth rate is 
linearly related to the amplitude of the parent mode. On the 
other hand, if 7 » |Ap|crcj, Re(s) ^ \Ap\'^a'^uj'^ /y, i.e., the 
growth rate is proportional to the square of the amplitude 
of the parent mode. The latter case is the one similar to the 
excitation mechanism we are discussing here. It should be 
noted that this parametric instability analysis gives a de- 
pendence of the growth rate in 7 which is not in agreement 
with the numerical results (considering 7 to be equivalent 
to l3), because the dependence of the spatial structure of 
the intermediate mode on the dissipation is not considered 
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in this simplistic analysis. Also, the parametric instability 
analysis suggests that the daughter modes gain energy from 
the parent mode, which is not what happens in the coupling 
mechanism we are considering, since here the differential ro- 
tation of the disc, and not the warp, is the ultimate source 
of energy for the r mode. Although the parametric instabil- 
ity analysis gives a dependence of the growth rate on the 
disturbance amplitude in agreement with our numerical re- 
sults, it is simplistic and does not consider all the details of 
the coupling mechanism. A parallel between the paramet- 
ric instability and our excitation mechanism is therefore not 
straightforward. 

In the strong dissipation regime, the intermediate mode 
dissipates completely, and does not influence the variation of 
the r mode growth rate with both the sound speed of the disc 
and the spin of the black hole. In this regime, the growth 
rate decreases with increasing Cg, as expected. The hotter 
the disc is, the wider the modes get, which means that if the 
sound speed is high, the modes are not as well trapped. More 
importantly, the shape of the warp changes when the sound 
speed changes since its wavelength (Aw) is proportional to 
Ca. The interaction relies on the use of the warp twice, there- 
fore we can argue that the growth rate is proportional to 
\dW/dr\'^ (since W represents the inclination and dW/dr 
the actual warp). Since |dVK/drp oc Wq /\w Wq /cg, i.e., 
for fixed sound speed the growth rate is proportional to the 
square of the inner warp amplitude (Fig.|4](b)) and for fixed 
Wo, the growth rate varies with (Fig. 11(c)). The smaU 
changes to the l/Cg law are justified by the fact that, for 
large sound speed, the decay rate due to the 'leakage' at 
rout is considerable. 

As for the variation of the growth rate with the spin 
of the black hole, an important conclusion is that, in fact, 
as argued in the beginning of this section, there is no mode 
excitation if the black hole is non-rotating. This was not ev- 
ident in Kato's simple estimates for the growth rates. The 
fact that the growth rate increases with a is also expected, 
not only because the waves are better confined for larger 
values of a, but mainly because the average warp ampli- 
tude in the trapped region is larger, for larger a. Also, the 
wavelength of the warp decreases as a increases, because 
the Lense-Thirring frequency (Olt = f2 — f2z « 2a/r^) in- 
creases. Therefore, for the same inner W, a larger dW/dr 
is achieved. In light of these very preliminary results, we 
could argue that HFQPOs would preferentially be detected 
in black hole candidates with large spin, if the interaction 
of the trapped r mode with a warped disc is the mechanism 
responsible for the excitation of the former. 



4.2 Growth rates in eccentric discs 

For the interaction between the r and intermediate modes 
in an eccentric disc, the variation of the growth rate with 
the dissipation factor is similar to the same variation for the 
interaction in a warped disc (Figs [7] (a) and|4](a)). Also, 
the change of the growth rate with the inner eccentricity 
(Fig. [7] (b)) is very similar to the variation of the growth 
rate with the warp tilt at the inner radius, i.e., there is 
a square dependence on the eccentricity amplitude at rin. 
This is expected as the eccentric mode plays, in the wave 
interaction described in this section, the role of the warp in 



the interactions described previously. Similarly, the variation 
with the sound speed (Fig. [7](c)) is also the expected one. 

The main difference between the interaction with the 
warp and the interaction with the eccentric mode is in the 
variation of the growth rate with the spin of the black hole 
(Fig. [7] (d)). The warp, being a n = 1 mode, has a variation 
with radius that strongly depends on the Lense-Thirring pre- 
cession frequency, and therefore, that strongly depends on a. 
On the other hand, the variation of the eccentricity ampli- 
tude with radius, given by equation (I37|l . is less dependent 
on the spin of the black hole. Therefore, the variation of a 
only causes variations of a factor of, at maximum, 2 in the 
growth rate obtained for the interaction within the eccen- 
tric disc. A very important difference is the fact that, in this 
interaction, a reasonable growth rate can be obtained in the 
case where a — 0. Therefore, in slowly rotating black holes, 
HFQPOs might be detected if the disc is eccentric, and if 
this excitation mechanism is responsible for the increase in 
the amplitude of oscillations. 



5 CONCLUSIONS 

In this paper we have described an excitation mechanism 
for trapped inertial modes, based on a non-linear coupling 
mechanism between these waves and global deformations 
(warping or eccentricity) in accretion discs. We have seen 
that the interaction of a trapped r mode with an intermedi- 
ate mode and a deformation in the disc results in growth of 
the trapped mode, if there is some process capable of mak- 
ing the intermediate mode dissipate in the disc. Dissipation 
is required so that this mode can remove rotational kinetic 
energy from the disc, which becomes available for the r mode 
to grow. Depending on the values of the sound speed, spin 
of the black hole and amplitude of the deformation at the 
inner radius, reasonable growth rates can be obtained for a 
warp or eccentricity of modest amplitude. In a warped disc, 
where the growth rate varies significantly with the spin of 
the compact object, growth rates as large as lo/10, where 
Lu is the oscillation frequency, can be obtained. If a = 0, 
no oscillations are excited in these discs. However, it may 
still possible to excite trapped modes in discs around non- 
rotating black holes if they are eccentric. 

The coupling process described here works as an exci- 
tation mechanism for trapped inertial waves, under a wide 
range of conditions, provided global deformations reach the 
inner disc region with non-negligible amplitude. The propa- 
gation of global modes, in a mor e realistic disc model, is the 
subject of a forthcoming paper (|Ferreira fc Ogilvi3l2008l ). 

In this paper we considered the excitation of trapped 
waves due to a non-linear coupling mechanism with global 
deformations, in a simple disc model. While this effect is re- 
sponsible for the growth of these modes, it has to compete 
with others that contribute to the damping of these waves. 
For example, since the conditions at the marginally stable 
orbit are unknown, it is possible for a 'leakage' of the trapped 
mode (similar to the one considered in the potential barrier 
analogy) through Vms to exist. This effect is to be considered 
in the future. Also, and more importantly, viscous dissipa- 
tion in the disc can cause damping of these modes. A simple 
estimate gives a damping rate of aQ. For small enough val- 
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ues of a and large enough warp or eccentricity, net growth 
can occur. 

Another point to be discussed is the applicability of 
our results to observed discs. We consider a very simple, 
isothermal disc model and wave perturbations for which 
7 = 1. In a more realistic disc, the vertical structure of 
the waves is changed while they propagate radially. The 
wave energy concentra tes either near the surface of the 
disc (jLubow fc Ogilvid IQQS,) or towards the disc mid-plane 
l|Korvcanskv fc Pringle 19951 ). which could potentially hin- 
der the propagation of intermediate modes away from the 
Lindblad resonance where they are excit ed. However, the 
proces s of 'wave channelling' mentioned bv lLubow fc Ogilvid 
l|l998f ) is only relevant at a distance from the resonance of 
~ ri^/m (rL being the radius of the Lindblad resonance), 
where the radial wavelength becomes comparable to the 
semithickness of the disc. Since the intermediate modes have 
m = 1, this effect is not important in the region where wave 
coupling occurs. The same is expected for cases in which the 
energy concentrates towards the disc mid-plane. The global 
deformation modes do not undergo significant wave chan- 
nelling because their wavelengths are always long compared 
to H. Therefore, we believe that our results, obtained in a 
simple disc model, are still qualitatively valid in more real- 
istic discs. 
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APPENDIX A: NON-LINEAR COUPLING EQUATIONS 
Al Wave coupling in a warped disc 

The following system of equations (cf. equations ([l}-((2|) describes the propagation of the r mode and n = intermediate 
mode, coupled by the warp, and needs to be solved for the growth rate to be determined: 



d 



d 



dt 



■KZ + "Rr - Sn^ik^ = --5— + Mr, 



dt 



dr 
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r d(f> 



mh'i + hi 



/r = {Mr, fB.4,, /rz) = -u! ■ VuW - Ww • Vul, 

/r;i = -wi ■ Vftw - ■ V/il, 

/i = (/ir, /:</., /iz) = -Wr ■ Vuw - "w • Vur, 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(A7) 

(A8) 
(A9) 
(AlO) 
(All) 



are the coupling terms, arising from non-linearities in the basic equations. 

Since we are interested in studying the axisymmetric r mode, the azimuthal mode number for the intermediate mode is 
1. Also, since the simplest possible r mode has one node in the vertical direction and the intermediate mode has n = 0, we 
use the following separation of variables: 



(«Rr,^R0,^R) = Re [(itRr(r-),'UR0(r-),/iR(r))Hei (^-^^ e~'"' 
Mrz = Re [■URz(r)Heo (^-^ e"'"*] , 

{uir,ui^,h'i) = Re |^(iiir(r),ui^(r),/ii(r))Heo (^-^) e"*^"'"*j , 
"^'iz = (2D mode), 

(uwri'Uw^i/iw) = Re \iuwrir),uw4,{r),hw{r)) ze''''^ , 
u'wz = Re l^uwz (r) c"^ j . 



(A12) 
(A13) 
(A14) 
(A15) 
(A16) 

(A17) 



This separation of variables results in having the coupling terms (|A8|I and (|A9|) proportional to Hei only, while the 
coupling terms (|A10|| - (|A11[) give rise to terms proportional to both Heo = 1 and He2 = z'^/H'^ — 1. Since we are interested 
in the terms that influence the mode with n = 0, we need to project these forcing terms on to Hcq. Also, the separation of 
variables results in coupling terms of the form 



Re{A)Re{B) = -Re{AB + AB* 



(A18) 



which means that the interaction of the m = r mode with the m = 1 warp results in two new modes, one with m = 1 and 
one with m = —1: 

r-mode (yl) oc e~'"* and warp (B) cx e"^ 

r-mode x warp oc + A_B* oc e"''""^* -f- e~"*~"^'. 

Since we are only interested in the action of this coupling on the intermediate mode with m = 1, because the one with 
m = — 1 does not have any Lindblad resonances (location where the wave should be excited) in the disc, these forcing terms 
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are projected on to e"* Similarly, the interaction of the intermediate mode with the warp gives rise to a m = 2 mode in 
addition to the ajdsymmetric r mode we are interested in: 

intermediate mode (A) oc e"*~''^* and warp (B) oc e"^ 

intermediate mode x warp oc AB + AB* oc e'2<^-i"* _|_ g~'"'_ 

Therefore, these forcing terms should be projected on to e~'"', which means that complex conjugates of warp quantities will 
appear in the equations. 

After separating variables, and projecting the forcing terms appropriately, the equations to be solved can be written as 



iLOUHr = 2iZuR0 ; — H + mi0— — H + ui^ H ■ — H - iu^^-—H (A19) 

dr 2 dr 2r r 2 dr ^ 2r 



. , (-.2 rr A{rUB..r) "ir dfe^ rr , • ^^W rr "Wr ^^1 „ . * /ll „ , . „„s 

-iLj/iR = il-^HuB.z — H + mi4,—-H —H - lUy^^ — H (A22) 

r dr 2 dr 2r 2 dr 2r 

r- , a\r-, , nn d/ll ItRr dUWr rr • '^Wr rr , rr '"Wz liWr dltRr rr 

- ILOUlr = -(l + pjilUlr + lilUl^ ■ ■ H - lUR^— H + UR0 H - Ur^ -— " " ■ H 

dr 2 dr 2r r 2H 2 dr 

-^iRz— (A23) 

, i/ll U-Rr duwth rr ■ ""Wt/) rr "Wr rr ^Wr dllR0 URr rr 

- lojtti^ = -(i + p)llui^ - —Ulr ; n - iu-R,4,— — H - ur0— — H ; — H - uw^-r—n 

2il r 2 dr 2r 2r 2 dr 2r 

-UR, - WW. ^ (A24) 

- luhi = + P)nhi Cs — — — — H - lUR^—H 7;-—^H - uy^-^ — - ur. — . (A25) 

r dr r 2 dr 2r 2 dr 2H 2 

This system is linear in the unknowns for the r and intermediate modes. The warp, which couples these modes together, is 
assumed to be known. 

For the interaction with the n = 2 intermediate mode, similar equations need to be solved. After separating variables 
and projecting forcing terms appropriately, the equations describing this interaction read 



DO dtiwr rr , • "Wr rr , n ^0 „ * duir . * Ulr ?tWz , ^Wr / . r,a\ 

iLUURr = 2i2uR0 ujr — ; -« + lujj, H + luis H - Uwr— ; — -n — iwwcft T^ujr H 7;—>J-iz (A26) 

dr dr r r dr ^ r H ^ 



K du-yy</) rr Uwr rr * '-"^It?) rr * '^ir „ "Wz , <(> / A o'7^ 

- IUUra, = -—rrURr - Itir — ; tl — Ui^ H - Uyfr— — H - Itwrf, T^^I^ H "Iz (A27) 

2S2 dr r dr r H 2 

-1C.UR. = -— - — - ^u^.^H - — ui. (A28) 

- ichn = f7^//nR. - £i^(!Z^ _ u.r^H + m,^^H - u%,r^H - iu^.^H - - (A29) 

r dr dr r dr r H 2 

r I «^0 I OO '^''I ^iRr dUWr rr • ""Wr „ , Uw^ „ UWr duRr / A nn\ 

- lujuir = - li + pjiluir + 2i2«i0 ■ ■ — H - iURd,— — H + URd, -H — H (A30) 

dr 2 dr 2r r 2 dr 

- lujui^ = -(i + f3)Qui^ - —Ulr -^H - mR4,—H - ur^—H - — —H - uw^ — H (A31) 

,. , 2hj URrduwz .uwz H duRz / a or,^ 

- ILUUlz = — (l + Pjiluiz t— U24, - —UWr ; (A32) 

U Z dr Zr Z dr 

■ U /• , a^O^, Cs d(riilr) 2i^I0 , o2 rr URr dhw „ ■ Wrr WWr d/lR / A OQ\ 

- luhi = -(1 + /3)nhi ^- — - - — ^ + Q^Huiz — H - iur^-—H — /f. (A33) 

r dr r 2 dr 2r 2 dr 

As before, the system is linear in the unknowns for the r and intermediate modes. It should be noted that although the same 
notation is used in the systems (|A19|| - (|A25|I and (|A26|) - (|A33|) to represent the intermediate mode quantities, they refer to 
two different modes: both with the same frequency and azimuthal mode number m — 1, but with different vertical mode 
number (n = in the first system and n = 2 in the second). 
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A2 Wave coupling in an eccentric disc 

After separating variables, and projecting the forcing terms appropriately, as done above for the interactions in a warped disc, 
the equations describing the coupling between the trapped r mode, eccentric disc and n = 1 intermediate mode read 

-iwmr^ = 2fiwR^ - — — — — ;r-uir + (iWBr + Me^) (A34) 

dr 2 dr 2 dr 2r 2r ^ 
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